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In several members of the ferro-pnictides, spin density wave (SDW) order coexists with supercon- 
ductivity over a range of dopings. In this letter we study the anomalous magnetic Zeeman response 
of this coexistence state and show that it can be used to confirm the extended s-wave gap structure 
as well as structure of superconducting (SC) gap in coexisting phase. On increasing the field, a 
strongly anisotropic reduction of SC gap is found. The anisotropy is directly connected to the gap 
structure of superconducting phase. The signature of this effect in quasiparticle interference mea- 
sured by STM, as well as heat transport in magnetic field is discussed. For the compounds with the 
nodal SC gap we show that the nodes are removed upon formation of SDW. Interestingly the size 
of the generated gap in the originally nodal areas is anisotropic in the position of the nodes over 
the Fermi surface in direct connection with the form of SC pairing. 

PACS numbers: 



The discovery of superconductivity at elevated tem- 
peratures in a new class of iron based materials (ferro- 
pnictides) [l] has revived interest in the underlying mech- 
anisms of high temperature superconductivity. Two key 
question are: what is the relation between antiferromag- 
netic order and superconductivity and what is the form of 
the superconducting pairing? These questions are funda- 
mentally connected to one another. A remarkable feature 
of the ferro-pnictides is that in many cases they exhibit 
coexistence of superconductivity and antiferromagntism, 
over a wide range of parameters. As discussed below, 
studying the coexistence phase can reveal important in- 
formation regarding both these questions. 

A popular theoretical model for superconducting pair- 
ing is the extended 5- wave {s±) pairing[2^-^, although 
other pairing symmetries were also discussed [TOl [11] . In 
its simplest form, s± pairing consists of nodeless singlet 
pairing but with different signs on different Fermi sur- 
face (FS) pockets. Experimental studies on the nature 
of pairing have not yet reached a unanimous conclusion. 
Hence new approaches to probing pairing are desirable. 

The parent compounds of many of the pnictides 
are antiferromagnetic metals [12], referred to here as 
the spin density wave (SDW) state. On doping, 
the magnetic order is reduced, and superconductivity 
emerges. Establishing the detailed phase diagram, and 
whether these two orders occur together, is an impor- 
tant question. In some members of pnictide family e.g. 
CeFeAsO, -cc-^ccQS] SDW and SC phases have no over- 
lap while in LaFeAsiO^F) [M", there are conflict- 
ing reports on SC and SDW phase coexistence. On the 
other hand coexistence of SDW and superconductivity 
in multiple pnictide materials have been reported. For 
example, in Bai-xKxFe2As2 and Ba{Fei-xCox) 2^82 
an extended region of coexistence with 0.2 < x < 0.4 
and 0.025 < x < 0.06 with a maximum superconducting 



transition temperature inside this region of ~ 28 and 
20 K respectively, are observed [161 [17]. Moreover, scan- 
ning probe measurements of the latter indicate that the 
two orders coexist in the same part of the sample [18 . 
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FIG. 1: a) FS in the first Brillouin zone, b) FS in the new 
Brillouin zone after SDW order forms. SC gap in red regions 
is robust against Zeeman. ko is the point where the gap closes 
first. c)SDW in transverse magnetic field 

Theoretical studies of SDW coexisting with super- 
conductivity has a long history (see [19 and references 
therein). For the pnictides, with extended 5- wave gap, 
the ordering wave vector Q = (tt, 0) [12^ folds the orig- 
inal Brillouin zone (figure |l(a)[ ) and parts of FS with 
opposite sign gap cross. It might seem that since SC 
gap changes sign on the new FS it has to have nodes. 
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But the sign change of the gap upon translation by Q 
(A^(q)A^(q + Q) < 0) protects the nodeless SC gap 
[I9l [20] . One simple way to understand this is that in 
the clean limit, it is possible to mathematically transform 
the s± SC with SDW problem, into an s wave SC with 
a charge density wave. Since the latter is an s-wave SC 
with a time reversal symmetric perturbation, its gap is 
protected by Anderson's theorem, and no nodes appear. 
This argument shows the possibility of coexistence of ex- 
tended s-wave SC gap and SDW but it is not exclusive. 
Indeed as we show simple 5-wave can also coexist with 
SDW despite the gap size reduces dramatically upon for- 
mation of SDW (also see [21 ) so coexisting phase alone 
is not indication of form of SC pairing. 

In this letter we study an unusual aspect of the re- 
sponse of the coexisting phase to a magnetic field which 
is direct evidence of extended 5- wave gap structure. The 
orbital effect of the magnetic field will be to produce vor- 
tices as in any type II superconductor. Here, we will focus 
instead on the Zeeman part of the coupling. In a layered 
superconductor which is appropriate for some of the pnic- 
tide materials, a field applied in the plane has a weaker 
orbital effect [22 . Moreover, the Zeeman coupling is ef- 
fectively enhanced by the presence of SDW order since 
those magnetic moments will cant along an applied field. 
The electron g factor is effectively increased 3 to 4 times 
over its bare value. The transverse ferromagnetic mo- 
ment that develops, can decrease the SC gap. We show 
that this gap reduction occurs in a highly nonuniform 
fashion. The gap remains robust in some regions of FS 
(red regions in fig. |l(b)[ ) up to very high fields. In these 
regions the pairing arises from the interplay of s± SC 
and SDW orders and effectively has a triplet character. 
So for the first time we show that coexisting phase can 
exclusively test the extended s-wave SC gap structure. 
Consequences for STM experiments as well as possible 
signatures of the Zeeman induced gap closing over parts 
of the FS in heat transport are analyzed. 

In some members of pnictides family, the SC gap has 
nodes and it is suggested that nodes are result of su- 
perposition of simple and extended s-wave pairing [23 . 
For the first time we study the coexisting phase of such 
pairing structure with SDW and show that upon forma- 
tion of SDW the nodes of this gap structure are removed 
with clear indication in multiple experimental probes e.g. 
heat capacity. The size of formed gaps are again strongly 
anisotropic in close connection with form of SC paring. 

Model: The Hamiltonian including coupling of the 
mean field SDW order with the conduction electrons is: 



Here = {i/j^^i/ji) and i/jg is the orbital spinor. For 
simplicity, in this letter we use a two band model for 
the pnictides [111 [25] which already captures the essential 
physics. So tjjj = {dxz-, dyz)g is the two component spinor 
corresponding to two d orbitals at the Fermi energy and 
the kinetic Hamiltonian is given bv[24]: 

h/c(k) = 2ti(cos kx — cos ky)X^ + 2(^2 — sin sin kyX^ 

+ [2(^2 + ^2) cos kx cos ky + 2t[{cos kx + cos ky) — jif]- \^ 

(2) 

where ti{t'i) and ^2(^2) represent nearest and next nearest 
neighbor hoping parameters, /i/ is the chemical potential 
and A is the Pauli matrix acting on the orbital space. 

Hspin corresponds to the spin interaction which leads 
to formation of Q = (tt, 0) ordering (e.g. Ji — J2 
model [26l [27]). Last term in ([T]) corresponds to coupling 
of ordered moments and the Zeeman field (H) with con- 
duction electrons (<j is the Pauli matrix acting on the 
Physical spin space). Saf and Sp are SDW moment 
and ferromagnetic moment generated by transverse field. 

Zeeman Coupling Strength: We now estimate the mag- 
nitude of the transverse Zeeman coupling (Ai?erro = 
\r]SF + f /i^Hl). Note, the SDW wih fiop into the plane 
perpendicular to the field, for sufficiently weak magnetic 
anisotropy, so the field is transverse to the SDW moment. 
To estimate the ferromagnetic moment Sp = xH, we use 
the transverse susceptibility in the non-SC SDW phase 
of the 1111 compounds: x 0-4 x 10~^ emu [27]. Since 
T] factor in ^ is not known, we can not estimate A 

Ferro 

directly. The ferromagnetic moment and SDW moment 
couple to the conduction electrons similarly so we can use 
the properties of SDW phase to estimate A Ferro- The 
magnitude of SDW moment is Saf ^ 0.36 /iB[12j and 
the SDW gap Asdw = \v^af\ ^ 0.08 eV^ is mea- 
sured using optical spectroscopy. The effective change in 
electron g factor is then 5g = j^r^^^^ - 5.5. The 
SC gap A5 ~ 4 meV is measured using ARPES[29]. In a 
transverse field of H 18T, Aperro ^ ^s- Iii the layered 
pnictides, a critical magnetic field of order 50T have been 
reported [22 so Zeeman coupling effect might indeed be 
as relevant as orbital effects in these fields. 

Coexistence Phase: To study the coexisting phase we 
use an extended Hilbert space and consider states at 
k and k + Q as two component pseudospin: = 
(^k, ^k-hq)- The Hamiltonian ([T]) with S^i^H^ and Sf\\x 
(fig. [T^ win be: 



Hmf = Yl ^kh(k)<l>k + Hspin 

k 

+ E (^^"^ ^AF + S^) 



(1) 



HMF(k) = ( ^'^^^ h,(k+Q) ' 

AsDW^^'T^O-^ + Aperro^^^^O-^ 



(3) 



r is the Pauli matrix acting on (q, q+Q)"^ space. Diag- 
onalizing h/e(k) we get the two bands. At each k point we 
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project in to the state closer to the Fermi energy which is 
a two component wave function '0s(k). Then h/g(k) will 
be replaced by the corresponding eigenvalue e(k)A^. Pro- 
jecting the Hamiltonian ([3| into the low energy orbital 
space we get: 



H^^(k) = [i?+(k)r° +i?-(k)r^] AV° 



^SDW 



A T a 



^Ferro 



\0 0x 
A T a 



(4) 



where E^{k) = £M±|]S±Q). when A 



Ferro 



the 

dispersion is e{k) = E^{k) ± v'£;-2(k) + A|^^. As- 
suming that the SDW ordering does not change the 
chemical potential, for the points on the FS we have 



e{k) = ^ ^+2(k) = E-^{k) + A2 



SDW 



The SC Hamiltonian acting on ^^k, ^-\^ is: 

H(k) = E+(k)rVV^ + £;-(k)rVV^ + A± (k)rV>^ 
+ A2^(k)rV>^ + AsDW /(k) rV>^ + AFerror'a^ fi' 

(5) 

where fi is the Pauli matrix acting on SC particle- 
hole space. A^^ is the the "simple" 5- wave (A^^(k) = 
A^^(k+Q)) and A^^ is the extended s-wave (A^^(k) = 
— A^^(k+Q)) part of pairing. We have traced over the 
orbital (tps) space leading to /(k) = ?/^J(k)?/^5(k + Q) 
in the SDW term which mixes the states at k and k + Q. 
This term is important in the perfect nesting limit and 
is ignored otherwise. 

First we consider Aperro = where the Hamil- 
tonian could be diagonalized analytically. Defining 
^2(k) = 2E+(k)2 + A±; + A^:, ^2(k) = ^2(k) _ 

Many features of coexisting phase could be understood 
from this dispersion: extended and "simple" s-wave 
pairing can both coexist with SDW although the gap 
for "simple" s-wave reduces greatly upon formation of 
SDW. If the SC gap has nodes as a result of presence of 
both "simple" and extended s-wave paring [23], upon for- 
mation of SDW the nodes are removed, i.e. magnetism 
enhances the SC properties! In general the gap formed 
is small [21 but is enhanced if nodes happen to occur 
close to the nesting regions. In the rest of paper A^^ = 
and use Age = Af^. 

Effect of Magnetic Field: Although it seems that all 
the symmetries of the Hamiltonian ([5| are broken by SC 
and SDW orders, there is a remaining symmetry imple- 
mented by the operator S = r^a^/i^ which commutes 
with the Hamiltonian in ([5|. S has four eigenstates with 
eigenvalue 1 and four with eigenvalue —1. We can reduce 
the size of the Hamiltonian in ([5| by projecting into the 
subspaces corresponding to different eigenvalues of S. 

For zero energy eigenvalue in the characterizing poly- 
nomial we can see that after the ferromagnetic moment 



is formed, at each k point the SC gap vanishes when: 



A2 A2 



E-ikY 



(6) 



Here one can readily see that the gap vanishes 
anisotropically, since E~(k) varies over FS, even though 



we consider magnitude of Asc to be uniform, 
the point where the SC gap vanishes first, the 
energy excitations dispersion is anisotropic: e{p 



At 
low 



E- 

the 



^0 



the 
the 



(ko)3 HPxl^PPy 

band curvatures, 
gap first closes (fig. |l(b) ) and p 
tion from ko. On the other hand 
when \AsDW^sc\ > E~{k.)'^ ferromagnetic moment 
can not close the gap! More specifically one can 
look at the point where ^~(k) = 0. Eigenvalues at 
this point could be calculated exactly. Defining i'^ = 



and (3 depend on 
the position where 

denotes devia- 
([6| shows that 



'^V^'sDW^'Ferro 



^SC^Ferro 



■^SDW^ energy is 



= + (A^^,,^ - A|^)2 + 4Asc/^SDw ± > 0. 

The eigenvalues are non-zero as long as both SDW and 
singlet superconductivity are present, regardless of the 
ferromagnetic moment. It is important to note that the 
gap is not SDW gap but it is indeed a SC gap (it vanishes 
when Asc = 0) that is robust against external magnetic 
field; as we will show it is indeed spin-triplet pairing gap. 

Two aspects of these results are particularly puzzling; 
since the singlet pairing changes sign on the FS, it seems 
that it should vanish at some points, but our result indi- 
cates that nodeless superconductivity coexist with SDW. 
The other feature is that the coupling with ferromagnetic 
moment affects the SC gap anisotropically and it can not 
destroy the gap in some regions. Below we show that the 
nature of the SC gap can explain these puzzling features. 

When Aperro = eigenvalues and eigenstates of the 
Hamiltonian for the points on the FS could be derived an- 
alytically. Interestingly two operators corresponding to 
the spin-triplet pairing also commute with the S. Op- 
erator that is important for us is TTripiet = t^g^ fi^ 
which anticommutes with the singlet operator Vsingiet = 
Asc^^o-y fj^y . The amplitude for any type of pairing could 
be calculated self-consistently using the wave functions 



l^n) 30 



^sc 



,E„<0 



{*„|rj,|«'„): 



As(k) CX ^(4'„|r<ji„g;et|*n) 

n 



(7) 



E- (k) (y^W+A|^ - E- (k)) 
^SDW - E-0<-) (x/i?-(k)2 + A|^^ - E-ik) 



'-^SDW 

Triplet 

n 

A-SDW 



= -Asc 



^Ai^^+E-ik^ 



As(k)2 +Ar(k)2 = A|c 



(8) 



(9) 
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As we expected A5'(k) vanishes where E (k) = (it 
is where singlet pairing changes sign) . Around this point 

A5(k) ^ Asc satisfies the expectation that 

singlet pairing should change sign between the regions 
coming from different FSs after zone folding. On the 
other hand when E~{'k) = 0, AT(k) ^ Asc so the pair- 
ing is triplet type. In the opposite limit E~ (k) ^ Asdw^ 
As{k) ^ Asc and AT(k) ^ A^c^^ so the pairing 
is mainly singlet. This even parity triplet pairing is ro- 
bust against coupling with the ferromagnetic moment i.e. 
regions with large triplet pairing remain gapped as the 
ferromagnetic moment forms. Eqn. (|9| also shows that 
singlet and triplet pairings together gap out all of the FS. 

A special limit "perfect nesting" (with t[ = 0)[24 . In 
this limit E~^{'k) = 0. It might seems that SDW gaps 
out all parts of the FS. But here /(k) plays an impor- 
tant role as it vanishes linearly at symmetry protected 
points on the FS[24 . With superconductivity a full gap 
opens, which closes at ko on increasing the field when 
\^Ferro\ = lA^d- The dispcrsiou then is 'semi-Dirac' 
like ^31j £^(ko + p) = y^^^f^p^^T^p^. We do not discuss 
this case further since it requires fine-tuning to reach. 

Experimental Consequences: So far we have proposed 
a theoretical picture to understand the coexistence of SC 
and SDW phase in ferro-pnictides. We showed trans- 
verse Zeeman field reduces the SC gap anisotropically. 
In the rest of this letter we will discuss the experimental 
signature of the effect discussed above. The usual experi- 
mental tool to map out the dispersion is ARPES which is 
not suitable in magnetic field. Another approach which 
by now is widely used to map the FS (e.g. in cuprate 
superconductors) is the quasiparticle interferences mea- 
surement using STM[32l [33]. STM measures the local 
density of state which is uniform for a normal clean metal. 
When sources of disorder such as impurities or crystal de- 
fects are present, elastic scattering mixes eigenstates that 
have different momentum but are located on the same 
contour of constant energy (STM bias voltage). When 
scattering mixes states ki and k2 , an interference pattern 
with wave vector q = ki — k2 appears in local density 
of states modulations which could be observed by STM 
as modulations of differential tunneling conductance [34 . 
The amplitude of the oscillation at momentum q is pro- 
portional to the joined density of state at momentum q 
i.e. n(k + q)n(q) where n(k) is the density of state at 
energy equal to the STM bias voltage and momentum k. 

The Zeeman field will generate the variation of the su- 
perconducting gap over the FS. The dispersion along the 
original FS is very shallow (at list of order Asdw^sc) 
compared to direction perpendicular to the FS. So the 
density of states at the tips of constant energy curves (fig. 

a)! is larger. The joint density of state for momentums 
connecting these regions (marked by arrows in figure 
is increased. These momentums vary continuously 
as STM bias voltage changes and will give the complete 



map of dispersion relation. Without external field, there 
is no gapless excitation. As external field is turned on 
the gap start to reduce anisotropically: vanishes in some 
regions but is not affected where E~ (k)'^ <C Asdw^sc- 
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FIG. 2: a) Parts of constant energy contours. The thick- 
ness of blue region corresponds to the density of states. 
b)Variation of SC gap over parts of the FS in transverse mag- 
netic field. Red and blue colors correspond the magnitude 
of triplet and singlet pairing respectively. The lines marked 
by arrows connect regions with large density of states at the 
energy marked by green curve. Fourier transform of density 
modulations have enhanced amplitudes at these momentums. 

Recent results on low temperature thermal conductiv- 
ity of pnictide superconductors [35, 36^ indicate the pres- 
ence of a full gap over the FS which is highly sensitive to 
the external magnetic field. The mechanism presented in 
this paper also leads to partial destruction of the SC gap 
in magnetic fields much smaller than the critical field. 
Based on our estimate using a uniform gap magnitude, 
this field is still much larger than the range where ex- 
periments have been performed. However, an anisotropic 
gap (as found in some calculations [6 ) could lead to much 
smaller onset fields where the ferromagnetic moment gen- 
erates nodes in the SC gap. When the gap first vanishes 
as external field increases, the low energy excitations have 
anisotropic dispersion e{p) cx a\px \ + which leads to 
density of states N{E) oc The signature of such 

density of states will be seen in temperature dependence 
of superfluid density Aps cx — a/T, as well as field de- 
pendence of heat conductivity Ahc (x (this could be 
understood as the Doppler shift due to superfluid flow 
around the vortex[37^). Note, in contrast a Dirac node 
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dispersion would have Aps cx — T and Atz oc . 

Conclusion: We considered the effect of Zeeman cou- 
phng in suppressing the extended 5-wave SC gap in the 
phase where superconductivity and SDW coexist and the 
transverse susceptibihty is enhanced. We showed that a 
highly anisotropic suppression of SC gap in Zeeman field 
is indeed directly related to the extended 5-wave struc- 
ture of SC gap. We acknowledge insightful discussions 
with Fa Wang and support from LBNL DOE-504108. 
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